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The latest XENON100 data severely constrains dark matter elastic scattering off nuclei, leading 
to impressive upper limits on the spin-independent cross-section. The main goal of this paper is 
to stress that the same data set has also an excellent spin- dependent sensitivity, which is of ut- 
most importance in probing dark matter models. We show in particular that the constraints set by 
XENON100 on the spin-dependent neutron cross-section are by far the best at present, whereas the 
corresponding spin-dependent proton limits lag behind other direct detection results. The effect of 
nuclear uncertainties on the structure functions of xenon isotopes is analysed in detail and found to 
lessen the robustness of the constraints, especially for spin-dependent proton couplings. Notwith- 
standing, the spin-dependent neutron prospects for XENON1T and DARWIN are very encouraging. 
We apply our constraints to well-motivated dark matter models and demonstrate that in both mass- 
degenerate scenarios and the minimal supersymmetric standard model the spin-dependent neutron 
limits can actually override the spin-independent limits. This opens the possibility of probing addi- 
tional unexplored regions of the dark matter parameter space with the next generation of ton-scale 
direct detection experiments. 



I. INTRODUCTION 

With hosts of on-going and forthcoming experiments, 
the direct search for weakly interacting massive parti- 
cles (WIMPs) is today a bustling field of research and 
one that has witnessed an impressive progress through- 
out the past years. The wealth of used target mate- 
rials and techniques means that we have now a rather 
constraining (but still inconclusive) picture of dark mat- 
ter scattering off nuclei. In particular, experiments us- 
ing proton- and/or neutron-odd target nuclei give useful 
limits on spin-dependent (SD) cross-sections, while de- 
tectors with high atomic number material chiefly probe 
spin-independent (SI) scattering - the former limits lag- 
ging significantly behind the latter. At the moment, the 
field stands at a stalemate with several collaborations in- 
cluding DAM A/LIBRA 1 , CoGeNT [2 and CRESST 
[3] hinting at a possible dark matter signal, a claim that 
has proven hard [U [5] to reconcile with the null results of 
XENON10/100 OE], CDMS [8] and others. The results 
of XENON10/100 have been of particular significance in 
challenging the existing hints since they rule out exten- 
sive regions of the dark matter parameter space. This 
extreme sensitivity is made possible by a combination of 
a very low background environment, a large exposure and 
the high atomic number of xenon isotopes. Recently [7J, 
the XENON100 collaboration has released its data corre- 
sponding to 225 live days and found essentially no events 
above the expected background. If interpreted in terms 
of SI dark matter elastic scattering, this non-observation 
leads to the world's best cross-section limits leaving all 
other experiments far behind (except perhaps in the low 
mass regime [5] [TU] or for isospin-violating dark matter 
candidates [H]). Additionally, as is well-known [51 112). 
detectors featuring xenon have also good sensitivity to 
SD scattering given the presence of the neutron-odd nu- 
clei 129 Xe and 131 Xe. It should be noted however that 



the structure functions of these isotopes do suffer from 
significant uncertainties [13HI9] that can affect the re- 
construction of dark matter parameters [20 . The aim 
of the present letter is precisely to point out that the 
current sensitivity level of XENON100 is already break- 
ing records in spin- dependent searches and to stress that 
upcoming ton-scale instruments are likely to deliver the 
strongest ever SD constraints in addition to their ac- 
claimed SI projected limits. In Section [TTJ we compute 
the limits set by the latest XENON100 data [7J on SD 
cross-sections and emphasise that these are highly com- 
petitive to SD-dedicated direct detection experiments. 
The effect of uncertainties on the structure functions of 
xenon isotopes is shown explicitly, and the prospects for 
XENON1T and DARWIN [2T] are presented as well. We 
illustrate in Section IIIII the usefulness of these SD con- 
straints by focussing on well-motivated particle physics 
models before concluding in Section [TV] 



II. SPIN-DEPENDENT SENSITIVITY 

In order to derive direct detection constraints we follow 
the standard computation of dark matter elastic scatter- 
ing rates in underground detectors [22] as implemented 
in our previous paper [25] . unless otherwise stated. We 
stick here to the so-called "standard halo model" [22] 
that features an isotropic Maxwell-Boltzmann velocity 
distribution of dark matter particles in our neighbour- 
hood, and use throughout a local dark matter density 
po = 0.3 GeV/cm 3 , a local circular velocity vq = 230 
km/s, a mean Earth velocity ve — 244 km/s and a local 
escape velocity v esc — 544 km/s. This set of values - 
chiefly based on the review by Lewin & Smith [22] and 
on |24j for v esc - is prone to astrophysical uncertainties 
[24] [25] . but has been widely used in the literature for 
the sake of comparison between different experimental 
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TABLE I: The spin expectation values (S p . n ) and structure functions Sij of the target nuclei 129 Xe and 131 Xe. 



results. In this framework the main uncertainty affecting 
the SD constraints from XENON100 regards the struc- 
ture functions Soo, Soi, Sn of the target nuclei 129 Xe 
and 131 Xe. We try to bracket this uncertainty and study 
its immediate effects by adopting the five parameteriza- 
tion specified in Tab. [j] "Bonn A" Q3], "Engel" |Hj and 
"Nijmegen II" Q~3j are rather standard nuclear models, 
while "Bonn CD" [TS] is more recent and "Menendez+" 
|19j refers to a very recent work. 

Unfortunately, we find some issues in quantifying nu- 
clear uncertainties in SD-proton scattering. Firstly, it 
appears that in the case of "Bonn CD" the combination 
of structure functions Soo + Sn + Soi (which is the rele- 
vant quantity for SD-proton scattering) almost vanishes 
for the energies of interest and is dominated by numer- 
ical errors. This corresponds to constraints on the SD- 
proton cross-section weaker by a factor ~ 1000 than those 
for the other parameterizations. Secondly, the "Menen- 
dez+" computations provide an error estimate for the 
structure functions Sqi and Sn which could be translated 
into an uncertainty on the SD limits. A naive combina- 
tion of these errors allows for a vanishing Soo + Sn + Sqi 
in the energy range of interest (i.e. no SD-proton limits), 
but these errors are probably correlated so that a sim- 
ple combination is not permissible. Since the appropri- 
ate prescription for including the errors is not available 
to us, we shall not consider them in the evaluation of 
the nuclear uncertainties. These two points indicate that 
the nuclear uncertainty in SD-proton scattering might 
be significantly larger than the canonical nuclear mod- 
els [T31 Q3] imply. Therefore, reliable limits on the SD- 
proton cross-section coming from xenon data require a 
better understanding of nuclear structure. None of these 
issues has a significant impact on the SD-neutron limits 
which depend on the combination Soo + Sn — Son 

Now, the key motivation for the present work is the lat- 
est XENON100 data set. In Ref. [7 the XENON100 col- 
laboration reports on the observation of two nuclear re- 
coil candidate events inside the WIMP signal region (en- 
compassing Er — 6.6 — 30.5 keV), while the background 
estimate amounts to 1.0 ± 0.2. These results correspond 
to a data taking period of 224.6 live days and an effec- 
tive exposure of 2323.7 kg. day. Applying the Feldman- 
Cousins procedure [25J with two observed events and 
a mean expected background of 1.0 events, we derive 
the 90% confidence level (CL) upper limit N R < 4.91 
for WIMP-induced nuclear recoils. Let us stress that 
more sophisticated procedures (e.g. a full profile likeli- 



hood analysis) would give somewhat stronger constraints, 
but our approach here is conservative and appropriate for 
this work. The above upper limit can be translated into 
separate constraints on the SD-proton and SD-neutron 
cross-sections ® , or into a combined constraint on both 
cross-sections using the approach of Ref. [27] (the stan- 
dard approach of Ref. [28 leads to weaker constraints 
than presented in Figs. [2] and [3] below). 

We present in Fig. [T] the 90% CL XENON100 upper 
limits on the SD-proton (left) and SD-neutron (right) 
cross-sections, together with the best experimental lim- 
its in the literature. The thick, short dashed, dotted, 
long dashed and dot-dashed red lines correspond to the 
five nuclear setups in Tab. |T] The nuclear uncertainties 
affecting XENON100 constraints can be rather large in- 
deed. As discussed in detail above, the SD-proton limits 
coming from xenon target experiments are extremely sen- 
sitive to nuclear computations (note that the "Bonn CD" 
<jp limit lies above the plotted range). Conservatively 
speaking, no solid upper limit on a^ D can be placed at 
the moment. In contrast, the nuclear nuisance on <7„ 
is relatively small and amounts to a factor ~ 2. Fig. [I] 
clearly shows that the latest XENON100 data are still 
not quite competitive in the SD-proton plane, lagging 
significantly behind SIMPLE and COUPP (that feature 
19 F), and especially IceCube. However, according to our 
analysis, XENON100 data beat the best published lim- 
its on the SD-neutron cross-section from XENON10 and 
ZEPLIN by approximately one order of magnitude, push- 
ing the upper limit down to cr„ < 5 x 10~ 40 cm 2 at 
m x = 50 GeV. 

Also shown in Fig. [T]are the prospects for XENON1T 
and a DARWIN-like xenon instrument, using in both 
cases the "Bonn A" setup. In the former case we sim- 
ply take a 60 times better sensitivity than XENON100 
(see [2T]), while in the latter instance we follow [35] and 
assume 2.00 ton.yr of effective xenon exposure for one 
background event and Er = 10 — 100 keV. This is a 
simplified approach which can be improved upon by the 
XENON1T and DARWIN collaborations themselves, but 
it is adequate for our purposes here. The prospects pre- 
sented in Fig. [I] are very encouraging. On the one hand, 
it seems feasible to constrain SD-neutron cross-sections 
of 10~ 42 — 10~ 41 cm 2 within the next decade. On the 
other hand, assuming "Bonn A", XENON IT will su- 
persede current SIMPLE and COUPP SD-proton limits, 
while DARWIN shall easily reach the present sensitivity 
of neutrino telescopes. However, current nuclear uncer- 
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FIG. 1: XENON100 90% CL exclusion limits for pure proton (left) and pure neutron (right) SD cross-sections. Shown are the 
results for five nuclear structure parameterizations: "Bonn A" 13 (thick red), "Engel" Q3] (short dashed red), "Nijmegen II" 
[T3] (dotted red), "Bonn CD" 18 (long dashed red) and "Menendez+" [19] (dot-dashed red). The prospects for XENON1T 
and a DARWIN-like experiment are indicated by the lower thick red curves. For a pure proton coupling (left) additional limits 
are available from IceCube [291 130| bb (empty blue circles), WW (filled blue squares), tt (filled blue circles), SIMPLE '12 |31) 
(black) and COUPP '12 (32] (magenta). For a pure neutron coupling (right) limits from CDMS '06 [33] (blue), XENON10 
'08 [5] (black) and ZEPLIN-III '12 [M] (magenta) are given for comparison. Notice that the "Bonn A" and "Engel" models 
lead to very similar results on a^ and thus the corresponding limits are indistinguishable, whereas the "Bonn CD" model 
yields a very weak Op D constraint above the plotted range. The light grey area represents the region of the mass-degenerate 
model as described in the text, while the dark grey region exhibits the measured relic density. Furthermore, we show a scan 
of the 11-parameter MSSM with the light blue dots; the green crosses indicate the points of the parameter space yielding the 
measured relic abundance. 



tainties are extremely large for SD-proton scattering on 
xenon. Only a better understanding of nuclear structure 
in the near future can eventually put XENON1T and 
DARWIN on the run for robust SD-proton limits. 

Recently, rather strong bounds on the SD-nucleon 
cross-section inferred from monojet searches at CMS and 
ATLAS 36-38 have also been presented, reaching the 
level of 10~ 40 cm 2 for axial-vector contact interactions. 
These constraints depend however on the effective op- 
erator type and can be considerably weakened for light 
mediators [39 . Similarly, neutrino bounds rely on the 
assumption of equilibration between WIMP capture and 
annihilation at the Sun. It is fair to mention at this 
point that neutrino telescopes [10], Large Hadron Col- 
lider (LHC) searches [55], direct detection experiments 
such as COUPP [H] and directional detectors [12] will 
all improve upon their present sensitivities along the next 
decade. In any case, the results in Fig. [T] give a strong 
motivation to take seriously the SD-neutron potential of 
ton-scale direct detection instruments. As we shall see 
in the next Section, the SD-neutron prospects of experi- 
ments such as XENON1T and DARWIN can easily over- 
shadow their SI prospects in the framework of different 



well- motivated dark matter models. 

Fig. |] displays the XENON100 90% CL limits on the 
SD effective couplings to neutrons and protons a„ iP for 
m x = 50 GeV, along with previous experimental results. 
The plot is particularly elucidative of the effect of nuclear 
uncertainties. It is nonetheless clear from Fig. [2] that 
XENON100 latest data shrink greatly the allowed region 
of parameter space, especially in combination with the 
results of COUPP and PICASSO. Finally, we show in 
Fig. [3] the combined SD constraints for m x — 50 GeV 
and a n /a p < 0. The choice of sign of a n /a p is motivated 
by the models discussed in the next Section. 



III. SPECIFIC CASE STUDIES 

We now illustrate the usefulness of the SD constraints 
from XENON100 in the framework of specific dark mat- 
ter models, starting with a simple mass-degenerate sce- 
nario [45Tl47| . This scenario consists of a Majorana dark 
matter particle x that couples to the right-handed up 
quarks via a coloured scalar i], C = —]x u bX\-, and is in- 
spired by a supersymmetric model with a squark that 
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FIG. 2: 90% CL exclusion region in the a„ - a p parameter 
space for a WIMP mass of 50 GeV. The line code is as in 
Fig.[l] We also show the limits from XENON10 '08 (blue), 
COUPP '08 [43] (black) and PICASSO '09 [44] (magenta). 
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FIG. 3: 90% CL exclusion region in the a s n D 
space for a WIMP mass of 50 GeV and a n /a p < 0. The line 
code is as in Fig. [T] Also shown is the theoretical expectation 
in the mass-degenerate model (grey) and in the 1 1-parameter 
MSSM (light blue). 



is nearly degenerate with the neutralino. In our analy- 
sis, we treat the coupling /, the mass m x and the ratio 
m v /m x > 1 as free parameters. In order to derive con- 
straints on the parameters of this model, it is first neces- 
sary to compute the WIMP-nucleon couplings from the 
WIMP-quark ones - for such we follow Ref. [23] and use 
the nuclear parameters therein. 

One interesting feature of this sort of models is that 
the mass degeneracy strongly enhances both SD and 
SI cross-sections. In fact, as shown in [33], even mod- 
est degeneracies lead to enhancements of several orders 
of magnitude. On the other hand, at sufficiently large 
WIMP masses and/or large degeneracies m ri /m x the to- 
tal event rate is dominated by the SD contribution, and 
thus the SD constraints from XENON100 become par- 
ticularly relevant. In order to show this point explic- 
itly, we perform a simple fixed-grid scan on the parame- 
ter space (f,m Xl m n /m x ), taking perturbative couplings 
in the range 10 -4 < / < 10 as well as m x > 40 
GeV, 10~ 2 < m n /m x — 1 < 10 2 and m, — m x > 1 
GeV. LHC constraints from jets and missing energy on 
coloured scalars that decay into first-generation quarks 
and a massive neutralino [48] apply only for relatively 
large mass splitting m v — m x > 10 2 GeV, while monojet 
and monophoton searches yield constraints only for rela- 
tively small splittings m v - m x < 10 - 20 GeV [4^1150] . 
Recently, it has been argued that a CMS razor analy- 
sis yields a constraint for mass splittings between 1 and 
100 GeV jSU] which, adapted to the model considered 
here, indicates a lower bound m n > 180 GeV. Note that 



LHC constraints [371 12H] on a^j D are based on contact in- 
teractions and cannot be applied directly, because of the 
light mediator rj. In particular, while the x _n ucleon scat- 
tering proceeds via the mediator in the s-channel, the cor- 
responding process at colliders would be a i-channel ex- 
change which is not significantly enhanced for ~ m x . 

Imposing the above collider constraints as well as the 
XENON100 limit on Si-only scattering, we get the light 
shaded region in Fig. [T] The dark shading indicates 
the region where, in addition, the observed dark matter 
abundance can be explained by thermal freeze-out when 
taking coannihilations into account [53J GO]. Clearly, a 
large portion of the parameter space escapes the usual 
XENON100 SI constraints but lies well above the corre- 
sponding SD-neutron upper limits. In other words, in a 
large portion of the parameter space the scattering rate 
in a xenon target experiment is dominated by the SD- 
neutron contribution. This stresses the importance of 
analysing XENON100 results in light of SD - and not 
only SI - scattering for mass-degenerate scenarios. It is 
also noteworthy that a chunk of the thermal region is 
already being excluded by XENON100 SD-neutron con- 
straints and that XENON1T and DARWIN will be able 
to probe a large portion of this region. The scan is also 
shown as a grey line in Fig. [3] for m x = 50 GeV. Note 
that it lies along a straight line because in these models 
a n /a p = Au(™V A « (p) - -0-52, or <J^ D /a^ D ~ 0.27. 

Let us now briefly comment on the prospects for 
supersymmetric neutralino dark matter. Within the 
constrained minimal supersymmetric standard model 
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(CMSSM), when taking LHC constraints into account 
and assuming a Higgs mass rrih — 125 GeV, most re- 
gions of the parameter space where the dark matter abun- 
dance can be explained by thermally produced neutrali- 
nos are highly constrained [52] (see also [551155) )■ In ad- 
dition, when applying constraints on SI scattering from 
XENON100, the SD cross-section gets pushed to values 
below 10 -42 cm 2 [53], which lie much below the future 
sensitivity of XENON1T or DARWIN (see Fig. [I). 

In order to check in how far this conclusion can be re- 
laxed when lifting some of the severe assumptions under- 
lying the CMSSM, we have performed an 11-parameter 
random scan of the MSSM using DarkSUSY [55] with pa- 
rameters shown in Tab. [Il] Apart from accelerator con- 
straints as implemented in DarkSUSY 5.0.5, including 
b — > S7 and the p-parameter, we require a Higgs mass 
in the range mt = 125.5 ± 1.5 GeV. Since a dedicated 
analysis of LHC constraints within the enlarged MSSM 
parameter space is beyond the scope of this work, we in- 
stead impose conservative bounds rriq 12 > 1.5 TeV and 
mq 3 > 500 GeV that are allowed within simplified mod- 
els [57] [55]. In addition, for each point we impose the 
XENON100 bound on Si-only scattering and conserva- 
tively require the annihilation cross-section to be below 
the most stringent Fermi-LAT limits from dwarf galaxies 
|59j . We have checked that within the range of param- 
eters, especially for vha and tan /3, recent constraints on 
BR(£? S — > [5U] are easily satisfied. However, al- 

lowing for larger values of tan /3 or a slightly wider Higgs 
mass range would not change our results significantly. 

The range of the SD cross-section that is compatible 
with all constraints is indicated by the light blue dots 
in Fig. [l] and by the light blue line in Fig. [3j Cross- 
sections as large as a^ D ~ 10 -39 cm 2 are reached for 
neutralino masses around 100 GeV - this is in line with 
the findings of Ref. [61]. The green crosses in Fig. [I] 
correspond to configurations yielding the measured relic 
abundance. The neutralino is in this case a mixed bino- 
higgsino state, characterised by a small //-term that is 
not much larger than Mi. Neutralino annihilation oc- 
curs mainly via Z and chargino exchange, while all pro- 
cesses involving squark exchange are suppressed due to 
the large squark masses. Note that the corresponding 
region within the CMSSM is disfavoured by LHC and 
Higgs mass constraints [52] . 

For the mixed bino-higgsino scenario, SI scattering is 
mediated mainly by Higgs exchange while SD scattering 
is mediated mainly by Z exchange, in the limit of heavy 
squarks [52]. Note that a n /a p ~ —0.8 for Z exchange, 
in good agreement with Fig. [3] In order to understand 
the relation between SD and SI scattering in this case, 
consider the effective Lagrangian 



(1) 



where the first (second) term describes SD (SI) scatter- 
ing, and i — u,d (in our numerical analysis we also take 
twist-2 and loop-induced couplings into account). In the 
decoupling/heavy squark limit, the coefficients for a neu- 



tralino x = z xi B + Z X? W + Z xZ H u 
by (we use the notation of [62] ) 



Z x ±Hd are given 
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The approximate expressions apply for Mi < \fi\ <C Mi 
and large tp = tan/3. Typically, both contributions are 
strongly correlated. Nevertheless, we would like to point 
out that in the 11-parameter scan we find a non- negligible 
portion of the parameter space where the SI cross-section 
is suppressed with respect to the SD one. Closer in- 
spection shows that this region corresponds to /z < 
and \n\/M 1 > 0(1). This behaviour is consistent with 
the approximate analytic expressions (see also [531 I51] L 
In particular, there exists a region in parameter space 
where a^ n < 10~ 46 cm 2 is below the Si-sensitivity of 
XENON1T, while the corresponding SD scattering cross- 
section on neutrons a^ D ~ 10 _40 cm 2 lies within reach of 
XENON1T. At colliders, an associated signal is a light 
chargino. Concerning monojet searches one naively ex- 
pects that the bounds inferred when assuming contact 
interactions get weakened by a factor of the order M§ /s, 
where s is the partonic centre of mass energy [39) . 



IV. CONCLUSION 

Direct dark matter detection is now entering a new sen- 
sitivity phase with different experiments - most promi- 
nently, XENON100 - pushing down cross-section upper 
limits to extremely low values. In this framework, it is 
important not to focus solely on SI scattering, but to ex- 
plore as well the SD potential of each data set. In fact, 
there are several dark matter models for which the SD 
contribution to the total event rate completely overshad- 
ows the SI one. As illustrated in this work, examples 
of such models are mass-degencrate scenarios and some 
parts of the parameter space of the MSSM. On that note, 
we have derived here the SD constraints set by the lat- 
est XENON100 data. While in SD-proton the constraints 
are prone to very significant nuclear uncertainties, in SD- 
neutron XENON100 overrides the best published limits 
by roughly one order of magnitude independently of nu- 
clear nuisances. Along the same lines, we find exciting 
prospects for XENON1T and DARWIN in constraining 
SD-neutron couplings. Interestingly, for a large class of 
dark matter models, the next generation of ton-scale di- 
rect detection instruments would detect first SD scatter- 
ing and only then the SI counterpart. 
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Parameter Mi M> f/7] m-L m^rt m^L m-R mj tan j3{Mz) rriA A~i 

Min [GeV] 10 80 80 1000 1000 500 500 80 5 500 -4000 

Max [GeV] 1 2000 2000 2000 3000 3000 3000 3000 3000 25 2000 4000 

TABLE II: Parameter ranges for the MSSM scan. The gluino mass parameter is fixed to M3 = 2 TeV, A Ut d,b = and 
mi = m/ R . All parameters except tan (3 and the pole mass are given at Q = 1 TeV. 
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